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1.  Introduction 


Fundamentals  of  thermodynamics  were  developed  during  two  centuries.  The  primary  notion  of 
thermodynamics  is  the  notion  of  the  equation  of  state  (EOS),  which  was  originally  associated 
with  something  later  coined  as  the  thennal  EOS.  Later,  it  was  changed  and  replaced  with  the 
concept  of  one  of  the  thermodynamic  potentials  referred  to  canonical  variables.  There  are  several 
reviews  describing  the  evolution  of  the  topic  and  its  current  state-of-the-art  (7-5). 

The  role  of  the  EOS  can  be  played  by  (1)  the  internal  energy  density  E  (V,  S) ,  referred  to  as  its 

canonical  variables — the  specific  volume  V  and  entropy  S  ;  (2)  the  free  energy  density 
F{V,T'),  referred  to  as  its  canonical  variables — specific  volume  V  and  absolute  temperature  T ; 

and  (3)  the  entropy  density  S(F,E)  ,  referred  to  its  canonical  variables — specific  volume  V  and 

the  internal  energy  E ,  etc. 

Finding  the  appropriate  EOS  is  very  important  for  any  discipline  using  continuum  description  of 
media.  Phenomenological  thermodynamics  suggests  some  canonical  experiments  to  determine 
the  required  thermodynamic  potential.  The  recommended  methodology  guarantees 
thermodynamic  consistency  of  the  model.  In  other  words,  this  methodology  guarantees  that  the 
model  obeys  the  First  and  Second  laws  of  thermodynamics. 

The  idealized  experiments,  recommended  by  theoretical  thermodynamics,  are  not  easy  or  even 
impossible  to  realize  practically  (e.g.,  when  dealing  with  planets  or  celestial  bodies).  Therefore, 
different  theoretical  and  experimental  techniques  for  establishing  the  EOS  have  been  suggested. 

A  combination  of  different  approaches  and  data  can  easily  violate  thennodynamic  consistency  of 
newly  suggested  models.  For  instance,  the  classical  models  of  the  ideal  and  van  der  Waals  gases 
are  thermodynamically  consistent  with  the  assumption  of  the  heat  capacity  Cv ,  depending  solely 

on  the  absolute  temperature,  and  are  incompatible  with  the  more  general  models  of  heat  capacity. 
As  it  was  demonstrated  in  Davis  (6),  the  classical  models  of  the  Mie-Griineisen  EOS  are 
thermodynamically  consistent  with  the  assumption  of  the  heat  capacity  Cv ,  depending  solely  on 

the  entropy  density  S  ,  and  are  incompatible  with  more  general  models  of  heat  capacity. 

Contrary  to  the  essentially  canonized  basic  principles  of  thennodynamics,  the  demand  for  novel 
EOS  always  remains  strong.  This  happens  because  of  the  pennanent  appearing  novel  substances 
or  because  of  the  necessity  of  working  within  extended  ranges  of  thermodynamic  parameters. 

Development  of  novel  EOS  is  unthinkable  without  further  progress  of  the  phenomenological 
thennodynamics  itself.  No  single  EOS  can  work  for  all  substances  and  all  ranges  of 
thennodynamic  parameters.  The  progress  and  revision  should  take  into  account  novel  types  of 
experiments,  not  just  those  idealized  experiments  considered  in  the  classical  textbooks. 
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Thermodynamics  and  theory  of  shocks  contain  strong  mathematical  cores  which  should  also  be 
further  developed.  The  models  appropriate  for  deep  mathematical  analysis  are  certain 
abstractions  of  the  practically  used  EOS.  By  sacrificing  some  details,  mathematical  analysis 
allows  a  much  deeper  analysis.  One  such  model  was  analyzed  in  Courant  and  Friedrichs  (7). 
We  suggest  its  generalization  in  the  next  section. 


2.  The  Suggested  Model 


The  generalized  Courant-Friedrichs  (gCF)  EOS  is  defined  by  the  following  relationship: 

E(V,S)  =  X(S)  Q(E)  +  A(E),  (1) 

where  V  is  the  specific  volume,  S  is  the  specific  entropy,  E  is  the  specific  internal  energy,  and 
X(S),Q(E),A(E)  are  three  functions.  The  functions  X(5'),Q(E),A(E)  should  be  specified 

based  on  available  experimental  data.  Thus,  we  arrive  at  the  EOS  based  on  three  functions  of 
one  variable  each — X(S'),  Q(E),  and  A(E). 


The  gCF  EOS  satisfies  the  following  differential  equation: 

JElnsMEl=0 

dVdS  8S 


(2) 


The  gCF  EOS  generalizes  the  Courant-Friedrichs  (CF)  EOS  with  the  separable  internal  energy: 


E(V,S)  =  X(S)  +  A(V), 


(3) 


which  is  specified  by  two  functions — X(S)  and  A (E)(7). 

The  CF  EOS  found  various  applications  in  modeling  shocks  in  condensed  liquid  media.  It  is  not 
only  useful  for  practical  needs  but  also  permits  deep  mathematical  analysis;  this  cannot  be 
overestimated.  As  it  was  obvious  from  the  very  beginning,  the  CF  EOS  has  limited  applicability. 
For  instance,  the  Griineisen  parameter  associated  with  this  model  vanishes  identically.  Because 
of  this,  we  think  it  makes  sense  to  analyze  the  properties  and  practical  applicability  of  the  gCF 
EOS  (equation  1). 

Using  equation  1,  we  arrive  at  the  following  formulas  of  the  absolute  temperature  T(V,S )  and 
pressure  P(V,S )  for  the  gCF  EOS: 

T{v,s)^n(v),  f(r.S).-x(S)«-«  ,4) 
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According  to  equation  4,  the  absolute  temperature  appears  to  be  the  product  of  the  functions  of 
5  andf. 


Also,  using  thermodynamic  identities,  we  arrive  at  the  following  formulas  of  the  heat  capacity  at 
constant  volume  Cv  (F,S)  and  the  Gruneisen  parameter T (V, S) : 

_*  Ata"!*!,  (5) 

Cr(V,S)  dS  dS 

and 

,  .  d  lnQ(F) 

nv^-v—±i.  <6) 

Equations  5  and  6  show  that  for  the  gCF  EOS,  the  heat  capacity  Cv  appears  to  be  the  function  of 
the  entropy  S  only,  whereas  the  Gruneisen  assumption  is  satisfied  automatically: 

cv  (v,s)  =  cv  (5),  r(v,s)  =  r(F) .  (7) 


It  was  demonstrated  in  reference  6  that  independence  of  T  of  5  and  independence  of  the  heat 
capacity  of  Cv  of  V  are  mutually  related  facts  because  of  the  universal  and  elegant 

thennodynamic  identity: 

a  r(r.s)  a  1 

es  v  evc,,(v,s)'  ' ' 

As  any  thermodynamic  function  of  state,  the  Gruneisen  function  (parameter)  can  be  expressed  in 
different  sets  of  thennodynamic  parameters.  It  can,  of  course,  look  different  depending  on  the 
choice  of  independent  thermodynamic  variables.  Among  the  many  presentations  of  the 
Gruneisen  parameter  are  the  following  two  (8): 

r (v,E)  =  vdF^V,E\  r {v,s)  =  -dXnT^V,S\  ( 9 ' 

V  ’  8E  V  ’  dlnV 

Per  the  terminology  of  reference  7,  equations  1  and  3  should  be  treated  as  the  complete  EOS, 
whereas  equations  5-7  can  be  treated  only  as  the  incomplete  EOS. 

In  the  case  of  constant  heat  capacity  Cv ,  we  get 


s-s * 

,  (10) 

and  the  gCF  EOS  can  be  presented  in  the  following  form: 
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(11) 


s-s * 

E(V,S)  =  Cve  Cv  Q(F)  +  A(F). 

Thus,  in  this  case,  the  gCF  has  only  two  functional  degrees  of  freedom. 


3.  Recovery  of  the  Functions  X(£),  Q(f) ,  and  A(f)  From  Experiment 


The  completeness  of  the  EOS  is  an  important  feature.  On  the  other  hand,  the  completeness  of 
the  EOS  in  no  way  guarantees  its  practical  usefulness.  It  only  guarantees  its  compatibility  with 
the  First  and  Second  Laws  of  Thermodynamics.  This  compatibility  is  important  but  not 
sufficient  for  ultimate  success.  The  EOS  can  be  useful  not  being  complete  and  can  be  complete 
not  being  useful. 

Real  usefulness  of  the  EOS  comes  from  comparing  it  with  experimental  data.  Fortunately,  the 
generalized  gCF  has  three  functional  degrees  of  freedom —  X  (S’) ,  Q(F) ,  and  A(F),  which 

matches  various  existing  experimental  data.  If  the  data  is  reliable  and  the  experiments  are 
chosen  wisely,  there  is  a  reasonable  chance  of  obtaining  useful  EOS. 

What  is  the  useful  set  of  experiments?  It  depends  on  the  intended  applications.  For  instance, 
measurements  of  the  heat  capacity  at  0  K ,  however  precise  they  are  and  however  useful  they  are 
for  the  low-temperature  physics,  can  be  quite  misleading  if  the  intended  application  deals  with 
modeling  a  supersonic  projectile  penetrating  through  a  thick  metallic  or  ceramic  wall  and  vice 
versa. 

For  determining  functions  X(S'),  Q  ( F ) ,  and  A(F) ,  each  set  of  experiments  requires  its  own 

formula  of  conversion.  Next,  we  provide  one  example  of  the  experiments,  giving  the  flavor  of 
the  typical  formulas  of  the  conversion.  These  fonnulas  are  basically  certain  integral  operators. 

Let  us  first  choose  the  reference  configuration  .  We  will  also  use  the  superscript  for 

all  quantities  relating  to  the  referenced  configuration.  Naturally,  in  order  to  recover  the  three 
functions,  we  need  three  sets  of  measurements.  We  will  consider  the  following  set  of  canonical 
experimental  measurements. 

Let  the  first  experiment  be  the  isochoric  heating  at  F  =  V* .  It  determines  the  heat  capacity 
Cv  (S)  at  constant  volume  as  a  function  of  the  entropy  density.  Let  the  second  experiment  be 

the  static  adiabatic  loading  corresponding  to  the  reference  value  S  -  S*  of  the  entropy  density. 
Let  T*  =  T*  (F)  and  F*  =  F*  ( F )  be  the  values  of  the  absolute  temperature  and  pressure  on  that 

very  adiabata. 
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Then,  the  measured  functions  Cr  ( 5'  ) ,  T*  =  T*  (F) ,  and  P*s  =  P*  (F)  determine  the  functions 
X(S'),  Q(F) ,  and  A  (E)  and  recover  the  corresponding  gCF  EOS.  The  recovered  equation 
reads  as  follows: 

s 

S  J  dljCy1  (ij)  V 

E(V,S)  =  T*s(V)\d£es *  -  j  dtuP* (m)  +  E*,  (12) 

s*  v* 

where  E*  is  the  internal  energy  density  in  the  reference  configuration. 

In  other  words,  the  fonnulas  of  conversion  for  this  set  of  experiments  read 

£ 

S  J*  drjCy1  (tj)  V 

X(s)  =  \d ^  ,  n(F)  =  T;(V),A(V)  =  jdmP;(ar).  (13) 

s*  v* 

In  the  case  of  constant  heat  capacity,  the  EOS  (equation  12)  can  be  rewritten  as  follows: 

S-S*  V 

E(V,S)  =  CX(V)eCv  -\dvjP;(m)  +  E*.  (14) 

V * 

When  the  Hugoniot  adiabata  data  (9,  10)  is  used  instead  of  the  data  relating  to  the  static  adiabata, 
the  inversion  formulas  (equation  13)  should  be  considerably  modified.  Additional  important 
aspects  should  be  considered  when  dealing  with  substances  undergoing  phase  transformations 
(see  references  9-11). 


4.  Discussion  and  Conclusion 


Following  the  footprints  of  Courant  and  Friedrichs  (7),  we  suggested  a  novel  complete  EOS 
called  the  gCF  EOS.  The  gCF  EOS  depends  on  three  functions  of  one  variable. 

For  the  case  of  two  available  experiments,  isochoric  and  adiabatic  loading,  we  established 
explicit  fonnulas  (equations  12  and  13),  which  converted  the  experimental  data  into  the  gCF 
EOS.  The  conversion  was  reduced  to  using  certain  integral  operations.  These  integral 
operations  are  stable  regarding  small  perturbations  of  experimental  data. 

The  thermodynamic  Griineisen  parameter  should  be  clearly  distinguished  from  the  Gruneisen 
assumption  and  the  Mie-Griineisen  EOS.  The  Gruneisen  parameter  for  any  two-parameter  media 
is  a  thermodynamic  function  of  state,  i.e.,  function  of  any  two  chosen  thermodynamic 
parameters.  It  is  a  meaningful  notion  applicable  to  any  two-parameter  media.  The  Gruneisen 
assumption  is  the  hypothetical  statement  applicable  to  a  limited  set  of  substances.  The 
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Griineisen  assumption  claims  that  the  Gruneisen  parameter  is  a  function  of  the  specific  volume 
only.  The  Mie-Griineisen  EOS  is  the  incomplete  EOS  based  on  the  Gruneisen  assumption. 

In  fact,  the  Mie-Griineisen  EOS  is  the  canonical  example  of  the  operational  EOS.  It  can  be 

presented  in  the  following  fonn: 

P(V,EpPr(V)  +  [E-Er(V)pf>-  ,  (15) 

where  the  functions  Py{V),  E  (V) ,  and  F (V)  are  functions  to  be  determined  from 

experimental  data.  It  is  the  example  of  the  incomplete  EOS  for  hydrocode.  Indeed,  as  any  EOS, 
it  is  a  function  of  two  variables — V  and  E  .  At  the  same  time,  equation  15  defines  an  operator 
transfonning  experimentally  measured  functions  P,(Vj  and  E  ( V )  into  an  incomplete  EOS. 

For  the  gCF  EOS,  the  Gruneisen  assumption  is  automatically  satisfied.  In  other  words,  the 
Gruneisen  parameter  appeared  to  be  the  function  of  the  volume  only,  whereas  the  heat  capacity 
appeared  to  be  the  function  of  the  entropy  density  only.  The  gCF  EOS  is,  in  fact,  the  most 
general  fonn  of  the  complete  EOS  that  is  compatible  with  the  Gruneisen  assumption. 

Other  examples  of  complete  EOS  will  be  given  in  the  forthcoming  papers  of  this  set  of  reports. 
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S  CROCKETT  MS  B221 
LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
MT  GREENFIELD  MSP915 
LOS  ALAMOS  NM  87545 


1  S ANDIA  NATIONAL  LAB  S 

S  SCHUMACHER  MS  0836 
ALBUQUERQUE  NM  87185 

1  S ANDIA  NATIONAL  LAB  S 

E  STRACK  MS  0378 
ALBUQUERQUE  NM  87185 

1  S ANDIA  NATIONAL  LABS 

T  VOGLER 

LIVERMORE  CA  94550 


LOS  ALAMOS  NATIONAL  LAB 
SR  GREENFIELD  MS  J565 
LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
R  MENIKOFF  MS  B214 
LOS  ALAMOS  NM  87545 


LAWRENCE  LIVERMORE  NATL  LABS 
L  BENEDICT 
M  A  BARRIOS 
A  TEWELDEBERHAN 
7000  EAST  AVE 
LIVERMORE  CA  94550 


LOS  ALAMOS  NATIONAL  LAB 
B  J  JENSEN  MS  B952 
LOS  ALAMOS  NM  87545 


NAVAL  SURFACE  WARFARE  CTR 
RESEARCH  AND  TECHLGY  DEPT 
G  T  SUTHERLAND 
INDIAN  HEAD  MD  20640 


LOS  ALAMOS  NATIONAL  LAB 

J  N  JOHNSON 

LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
B  PLOHR  MS  B213 
LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
L  SMILOWITZ  MS  P915 
LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
D  PRESTON 

LOS  ALAMOS  NM  87545 

S ANDIA  NATIONAL  LABS 
J  H  CARPENTER  MS  1322 
ALBUQUERQUE  NM  87185 

S ANDIA  NATIONAL  LABS 
M  DESJARLAIS  MS  1189 
TR  MATTSON  MS  1189 
ALBUQUERQUE  NM  87185 

S ANDIA  NATIONAL  LABS 
RJ  MAGYAR  MS  1322 
ALBUQUERQUE  NM  87185 


1  APPLIED  RESEARCH  ASSOC 
D GRADY 

4300  SAN  MATEO  BLVD  A-220 
ALBUQUERQUE  NM  871 10 

1  US  AIR  FORCE  RESEARCH  LAB 
MUNITIONS  DIRCTRT 

Y  HORIE 

EGLIN  AFB  FL  32542 

2  THE  JOHNS  HOPKINS  UNIV 
DEPT  OF  MECH  ENGRG 

K  T  RAMESH 
T  W  WRIGHT 
232  LATROBE  HALL 
3400  N  CHARLES  ST 
BALTIMORE  MD  21218 

2  UNIVERSITY  AT  STONY  BROOK 
J  GLIMM 
R  SAMULYAK 

STONY  BROOK  NY  11794-3600 

2  UNIVERSITY  OF  SOUTH  FLORIDA 
DEPT  OF  PHYSICS 
I  OLEINIK 

V  ZHAKHOVSKY 
TAMPA  FL  33620 
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COPIES  ORGANIZATION 


1  UNIV  OF  UTAH 

DEPT  OF  MECH  ENGRG 
R  BRANNON 

50  S  CENTRAL  CAMPUS  DR 
SALT  LAKE  CITY  UT  841 12 

1  UNIV  OF  CALIFORNIA 

DEPT  OF  MECH  &  ARSPC  ENGRG 
V  NESTERENKO 
9500  GILMAN  DR 
LA  JOLLA  CA  92093-0411 

1  UNIV  OF  MD 

DEPT  OF  MECHANICAL  ENGRG 
R  ARMSTRONG 
COLLEGE  PARK  MD  20742 

1  UNIV  OF  COLORADO  AT  BOULDER 
DEPT  OF  CIVIL  ENGRG  & 
ARCHITECTURAL  ENGRG 
R  A  REGUEIRO 
428  UCB  ECOT  441 
1 1 1 1  ENGRG  DR 
BOULDER  CO  80309-0428 

1  ROBERT  R  MCCORMICK  SCHOOL 
OF  ENGRG  &  APPLIED  SCI 
H  ESPINOSA 
633  CLARK  ST 
EVANSTON  IL  60208 

1  UNIV  OF  MISSOURI 
DEPT  OF  CHEMISTRY 
T  SEWELL 

COLUMBIA  MO  652 1 1  -7600 
1  UNIV  OF  MD 

DEPT  OF  MATERIALS  SCI  &  ENGRG 
M  KUKLJA 

COLLEGE  PARK  MD  20742 

1  STANFORD  UNIV 

DEPT  OF  AERONAUTICS  & 
ASTRONAUTICS 
R  CHRISTENSEN 
DURAND  BLDG  496 
LOMITA  MALL 
STANFORD  CA  947305-4035 

1  US  ARMY  RESEARCH  OFFICE 
RDRL  RO 
D  SCUTRUD 
BLDG  4300 
DURHAM  NC  27703 


1  US  ARMY  RESEARCH  OFFICE 
RDRL  ROI 
B  WEST 
BLDG  4300 
DURHAM  NC  27703 

1  US  ARMY  RESEARCH  OFFICE 
RDRL  ROI  M 

J  MYERS 
BLDG  4300 
DURHAM  NC  27703 

2  US  ARMY  RESEARCH  OFFICE 
RDRL  ROE  M 

J  PRATER 
D  STEPP 
BLDG  4300 
DURHAM  NC  27703 

1  US  ARMY  RESEARCH  OFFICE 
RDRL  ROE  N 
R  ANTHINIEN 
BLDG  4300 
DURHAM  NC  27703 

3  US  ARMY  RESEARCH  OFFICE 
J  LAVERY 

D  MANN 
S  MATHAUDHU 
BLDG  4300 
DURHAM  NC  27703 

ABERDEEN  PROVING  GROUND 

72  DIR  USARL 
RDRL  WM 
B  FORCH 
J  MCCAULEY 
SKARNA 
P  PLOSTINS 
RDRL  WMP  B 
R  BECKER 
S  BILYK 
D  CASEM 
J  CLAYTON 
D  DANDEKAR 
J  FITZPATRICK 
M  GRINFELD  (20  CPS) 

C  HOPPEL 
R  KRAFT 
B  LEAVY 
D  POWELL 
M  RAFTENBERG 
M  SCHEIDLER 
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T  WEERASOORIY A 
C  WILLIAMS 
RDRL  WMP  C 
T  BJERKE 
S  SEGLETES 
W  WALTERS 
RDRL  WML  B 
I BATYREV 
J  BRENNAN 
S  IZVEKOV 
B  RICE 
RDRL  CIH  C 
P CHUNG 
J  CAZAMIAS 
J  KNAP 
RDRL  WMM  G 
I  ANDZELM 
F  BEYER 
J  LENHART 
B  R1NDERSPACHER 
RDRL  WMP  E 
W  GOOCH 
B  LOVE 

C  KRAUTHAUSER 
RDRL  WMM  E 
J  ADAMS 
L  KECSKES 
J  LASALVIA 
P  PATEL 
J  SWAB 
M  WILL-COLE 
RDRL  WMM  B 
B  CHEESEMAN 
C  FOUNTZOULAS 
G  GAZONAS 
RDRL  SER  L 
W  NOTHWANG 
RDRL  WML  H 
M  FERMEN-COKER 
RDRL  WMP 
B  BURNS 
S  SCHOENFELD 
RDRL  WMP  D 
H  MEYER 
RDRL  WMP  G 
R  BANTON 
S  KUKUCK 
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COPIES  ORGANIZATION 


1  CAVENDISH  LAB 

MOTT  BLDG  RM316B 
JJ  THOMSON  AVE 
M  CHAUDHRI 
CAMBRIDGE  CB3  OHE 
UNITED  KINGDOM 

1  LOUGHBOROUGH  UNIV 

WOLFSON  SCHOOL  OF  MECH  & 
MANUFACTURING  ENGRG 
LOUGHBOROUGH 
V  SILBERSCHMIDT 
LE11  3TU 

UNITED  KINGDOM 

1  UNIVERSITY  OF  LIVERPOOL 
DEPT  OF  ENGRG 
BROWNLOW  HILL 
LIVERPOOL 
OUYANG 
L69  3GH 

UNITED  KINGDOM 

4  AWE  ALDERMARSTON  READING 
BERKSHIRE  RG7  4PR 
N  BOURNE 
J  C  F  MILLETT 
G  A  COX 
C  M  ROBINSON 
UNITED  KINGDOM 

1  DEFENCE  SCI  &  TECHLGY 
ORGANISATION 
WEAPONS  SYSTEMS  DIVISION 
A  RESNYANSKY 
EDINBURGH  SA5111 
AUSTRALIA 

1  CSIRO  EXPLORATION  &  MINING 
PO  BOX  883 
KENMORE  QLD  4069 
F  D  STACEY 
AUSTRALIA 

1  MOSCOW  STATE  UNIV 

INSTITUTE  OF  MECHANICS 
MICHURIN SKI(E  PR  1 
S  S  GRIGORYAN 
MOSCOW  117192 
RUSSIA 


1  STEKLOV  MATHEMATICAL 
INSTITUTE 
GUBKINA  AV  8 
A  G  KULIKOV SKII 
MOSCOW  117966 
RUSSIA 

1  LANDAU  INSTITUTE  FOR 
THEORETICAL  PHYSICS 
N  INOGAMOV 

RAS  CHERNOGOLOVKA  142432 
RUSSIA 

5  JOINT  INSTITUTE  FOR  HIGH 
TEMPERATURE  PHYSICS 
RUSSIAN  ACADEMY  OF  SCIENCES 

V  FORTOV 

V  GRYAZNOV 
GKANEL 

M  LOMONOSOV 

S  RAZORENOV 

JIHT  RAS  MOSCOW  125412 

RUSSIA 

1  INSTITUTE  OF  CONTINUOUS  MEDIA 
MECHANICS 
URAL  BRANCH  OF  RAS 
O  B  NAIM  ARK 
1  ACAD  KOROLEV  STR 
PERM  614013 
RUSSIA 

1  BEN  GURION  UNIVERSITY 

DEPARTMENT  OF  MECHANICAL 
ENGRG 

E  B  ZARETSKY 
PO  BOX  653 
BEER  SHEVA  84105 
ISRAEL 

1  RAFAEL 
PO  BOX  2250 

V  PARTOM 
HAIFA  3 1021 
ISRAEL 
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